Precessional motion of a vortex in a finite-temperature Bose-Einstein condensate 
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We study the precessing motion of a vortex in a Bose-Einstein condensate of atomic gases. In 
addition to the former zero-temperature studies, finite temperature systems are treated within the 
Popov and semiclassical approximations. Precessing vortices are discussed utilizing the rotating 
frame of reference. The relationship between the sign of the lowest excitation energy and the 
direction of precession is discussed in detail. 
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I. INTRODUCTION 

The Gross-Pitaevskii (GP) approximation is often 
used to analyze the Bose-Einstein condensates (BEC) 
of atomic gases fj] • The GP approximation only treats 
the condensate fraction, leaving out the normal com- 
ponent. Therefore, the GP equation treats systems at 
zero temperature. The Bogoliubov equations [l| describe 
the excitations supported by the condensate in the zero- 
temperature limit. The Bogoliubov excitation spectra 
agree with experimental results, such as the collective 
oscillation modes (including the Tkachenko waves) 13, Q 
of a vortex lattice and those on Bragg spectroscopy |J| of 
elongated condensates. Since the physical systems exist 
at finite temperatures, attempts to incorporate effects of 
temperature into the various approximations have been 
made [a, |(j . Among these, the Hartree-Fock-Bogoliubov- 
Popov (Popov) approximation (Sec. V. F in Ref. 0) is 
one of the most common finite-temperature approxima- 
tions. 

The Popov approximation treats the density of the 
normal gas component as a mean field. The condensate 
is described with the GP equation, extended to include 
the mean field potential. The excitation spectra and the 
wavefunctions of the normal gas are given by eigenequa- 
tions which are similar to the Bogoliubov equations. The 
excitation spectra and the wavefunctions constitute the 
density of the normal gas. 

Within the Bogoliubov (T = 0) approximation, the oc- 
currence and the disappearance of a vortex are indicated 
by the existence of an excitation with negative excita- 
tion energy (in other words, the anomalous excitation). 
There also exist relations between the direction of the 
precessional motion of the vortex and the sign of the ex- 
citation energies within zero-temperature theories. This 
paper aims to answer the following question: Does this 
relation stay valid in the Popov (T > 0) approximation? 

Besides the problem of the direction of precession, the 
precession frequency within the Popov approximation 
has not been investigated in detail. The precession fre- 



quency at zero temperature has been discussed within the 
Thomas-Fermi (TF) approximation using several meth- 
ods (see Sec. 5.1 of Ref. @ and Ref. 0). There also 
exists an analysis within the classical-field approxima- 
tion |lC| in finite-temperature systems. 

The condensate is assumed to be trapped with a rota- 
tionally symmetric trap. Therefore, a precessing vortex 
is an off-centered vortex. Within zero-temperature the- 
ories, a static system with an off-centered vortex ,11] in 
the rotating frame is equivalent to a system in the pres- 
ence of a precessing vortex. If the displacement Ar of 
the vortex from the trap center is small, the excitation 
energy of the anomalous mode 01 > or the core-localized 
mode depends linearly on the rotation fre quen cy. The 
variational Lagrangian analysis used in Ref. [12] and Bo- 
goliubov theory yield the same dependence for a centered 
vortex. 

In a previous study , we discussed a precessing vor- 
tex state as a deviation from the axisymmetric configura- 
tion. In this work we aim to treat the precessing vortex 
in a 2D geometry directly within the Popov approxima- 
tion in order that the relation between the precession and 
the excitations are directly taken into account from the 
outset. 



II. APPROXIMATIONS 

The condensate is treated with a nonlinear Schrodinger 
equation (NLSE) within the Popov approximation. The 
thermal atoms are described using the Popov equations 
which are eigenequations. Because of computational 
complexity, we use the Popov equations only for exci- 
tations below a certain cutoff energy i? cu t. The exci- 
tation energies e > E cut are taken into account within 
the semiclassical approximation 0, 0, The den- 
sity of the condensate is n (r) = \<j)(r)\ 2 , where 4>{r) 
is the condensate wavefunction. The normal particle 
density coming from the excitations below (above) E cxli 
is 7ii (^2). Therefore, the particle number density is 
n(r) = no(r) + ni(r) + ri2(r). The condensate is de- 
scribed with the NLSE 
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{-CV 2 +V-[j,+g(n a + 2n 1 + 2n 2 )-u i0 frxp}(f> = (1) 
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where C — Ti 2 /(2m) and g = AtkU a/m. The mass of 
a Na atom m = 38.17 x 10 - 27 kg, and the scattering 
length a = 2.75 nm for Na atoms are employed. We use 
the cutoff energy E cut — lOftwtr E3- Angular velocity 
of rotation is Lu mt and the rotation axis is parallel with 
the z axis (u; ro t = e z 07 ro t). Excitations below the cutoff 
-E cu t are eigenstates of the Popov equations 

{-CV 2 + V - [J,+ 2g(n + n x + n 2 ) 

-w rot • r x p}u q - g4> 2 v q = e q u q , (2a) 

-{-CV 2 + V - n + 2g(n + m + n 2 ) 

+uj rot -r xp}v q + gcj)* 2 u q = e q v q . (2b) 

These reduce to the Bogoliubov equations if we neglect 
ni and n%. The wavefunctions u q and v q obey the nor- 
malization condition 



(h 



dr = 1. 



(3) 



The density n\ is a weighted sum of the wavefunctions u 
and v: 



" 1 1 r 1 = < E (KI 2 + KI 2 )/( £ «) + KI 2 K 4 ) 

_ 9 (e,<-E cut ) 
1 



/(e) = 



#bt) - r 



(5) 



The higher-energy range e > Ecnt is described within the 
semiclassical approximation 0, ITHl ITfil ] which neglects 
the derivatives of the amplitudes of the wavefunctions u 
and v and the second derivatives of their phases. We also 
neglect the phase of the condensate wavefunction <j> here. 
Then the Popov equations reduce into algebraic form 
(Ref. [16| , Eqs. (5)). The expression for n 2 , in analogy 
with Eq. is: 



^ = /f (/(*> + 5) " 5 W'- 1 ^ 

where the Hartree-Fock (HF) energy 



(6) 



£hf(t",p) = ^ 1- V - M + 2 .g(n + n : + n 2 ) (7) 

2m 



and energies 



e HF( r ;P) - .9 2 «o, 



e(r,p) = 
e(r,p) = i(r,p)-v Tot -rxp 



(8) 
(9) 



are functions of r and p. The noncondensate densities ni 
and ri2 are determined from Eqs. Q and ©. They are 
treated as mean field potentials throughout the above 
equations. Thus the numerical procedure needs to be 
selfconsistent such that it is repeated until the solution 
reaches convergence in which all the equations are simul- 
taneously satisfied. 



The angular momenta of the condensate (j> and of the 
wavefunctions u, v are 

y . f 6*{r X p)cbdr , x 

?i J no(r)ar 

At \ Ju*(rxp)u q dr 
Aug) = e z ^-^^ , (11) 



fiU a 



A(v q ) = e z 



Jv*(r x p)w,jdr 



(12) 



where U q = f \u q \ 2 dr and V q = J \v q \ 2 dr. Within the 
Bogoliubov theory (T = 0), the excitation energies de- 
pend linearly on the angular velocity uj mt as follows: 

£ = eiab - TiLo rot q e , (13) 
qe ee Rc [{A{u q ) - A{0)}U g + {A(v q ) + A(&}V q ] 
/{U q + V q ). (14) 

The value qe is useful for characterizing 0, Q an excita- 
tion also for finite-temperature systems. But as for the u> 
dependence of e, there are deviations from Eq. (|13f) due to 
changes of density in the normal component n\(r) +n2(»") 
for finite-temperature systems. Figure indicates these 
deviations. Changes in the excitation energies modify 
the normal-component density which in turn affects the 
whole system, including the e's themselves. Therefore, 
the spectra for each value of w ro t needs to be calculated 
individually. 



III. OFF-CENTERED VORTEX 

The sign of the excitation energy of the core-localized 
excitation and the direction of the precessional motion 
are related flU ] at zero temperature. The predicted pre- 
cession frequency fits well with results of the experi- 
ments 18] . A finite-temperature extension of the Bo- 
goliubov equation, the Popov approximation, shows that 
the sign of the core-localized excitation becomes posi- 
tive |19l l20j . If the direction of the precessing motion 
and the sign of the excitation energy correspond to each 
other, the precessing motion of the vortex must also be 
inverted as the lowest excitation energy raises from neg- 
ative to positive values. 

We extend the 2D treatment of the Bogoliubov 
equations to finite-temperature Popov equations. It 
makes possible to treat the slightly off-centered vortices 
directly within the Popov approximation. 

Assume a BEC system has a vortex line perpendicular 
to the z axis. Particles are confined with a harmonic trap 
along the x and y axes 



V(x,y) = 



-(* 2 + y 2 ) 



(15) 



with w tr = 2ir x 200 Hz. The system has finite z thick- 
ness and it is uniform along the z axis. Periodic bound- 
ary conditions along z are used. Hence this system is 
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FIG. 1: Dependence of e' (lowest e) on w ro t in an axisymmet- 
ric system (Ar = 0). The solid lines indicate the e at T = 
(Bogoliubov approximation) and T/T c = 0.1,0.3, and 0.5 
(Popov approximation). The dotted lines represent Eq. 113ft . 
The angular momentum qg — —1. The dotted line and the 
solid line overlap at T = 0. The two lines are separated for 
T > 0, which means that "equality" between the rotation 
velocity w ro t and the excitation energy in Eq. 113H is lost at 
finite temperatures. 



not a two-dimensional BEC 21, 22], but rather a three- 
dimensional BEC having restricted geometry. We treat a 
system having 10 5 atoms within a z thickness of 10 /im. 
The Thomas-Fermi (TF) radius Rtf = 6.793 (im is used 
as the scale of length. 

Equations JTJ - J3J are repeatedly solved until conver- 
gence into a self-consistent solution. While the vortexfree 
(A{4>) = 0) and centered-vortex {A{<p) — 1) configura- 
tions are most likely, there also exists a solution with an 
off-centered vortex (0 < A((f>) < 1) in a narrow window of 
u>. Figures [21 a-c) represent such a typical system. Here 
the angular momentum is A((f>) = 0.863. The particle 
number of the condensate is 48% of the total particle 
number. The noncondensate density has a characteristic 
peak at the core of the vortex, like those in the axisym- 
metric studies 0, |2(J • 

The displacement of the vortex core Ar is unrestricted 
in the numerical processes, unlike for the axisymmetric 
situations (Ar = 0). Therefore, the displacement Ar 
depends on temperature T and the rotation frequency 
Lu mt as presented in Fig.[3ja). 

Figure[3f a) plots rotation frequencies at which the sys- 
tem is static. Let us denote the displacement and the 
rotation frequency at the static point as Ar' and oj' Iot . 
When Ar < Ar', w rot = tj' Iot , and T = 0, the vortex 
has an instability and it tends to move inward. When 
Ar > Ar', the direction is outward. As w ro t increases 
in Fig. Ota), the system has a wider range of displace- 
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FIG. 2: Density profiles of a system at the temperature 
T = 0.5 T c , and angular momentum of the condensate A(4>) = 
0.863 in rotating coordinates with u} ro t = 0.300 Ufa. (a) Par- 
ticle density no(x,y) of the condensate, (b) Particle density 
ni(x,y) +ri2(x,y) of the noncondensate, (c) Excitation spec- 
tra e q vs. qg. The solid line shows 1 — Lu ro tqg / uitr ■ Its slope in- 
dicates the rotation velocity ui ro t ■ The line touches the dipole 
modes at q$ = ±1. 



ments Ar having an inward instability. This instability 
brings the vortex to the center of system and makes the 
vortex state more sustainable. At the finite temperature 
T = 0.1T C , the range of Ar having an inward instability 
is almost the same. But the positive value of the lowest 
excitation energy e' in Fig. B^b) requires that the dis- 
placement Ar is stable and the vortex tends to remain 
at Ar = Ar'. 

The lowest excitation energies e' in Fig.|3Jb) differ sig- 
nificantly between the T — system and those of the 
T = 0.1 T c system. It is small and negative (0 > e' > 
— 0.004?kJtr) for the T = system, while it is positive 
in the T = 0.1 T c system. The corresponding rotation 
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frequencies w rot in Fig. OJa) differ only a little. This 
indicates that the direction of the precessional motion 
which is represented by the sign of w rot and the sign of 
the lowest core excitation energy are not related in the 
T > systems. This point is discussed further in the 
next section. 

Figures |2fa-b) compare the spectral densities 

.91 (i) = Jdr (M 2 + KI 2 )/(^) + N 2 > 

(16) 

^) = /^{^(/(,) + i)-i}^- £ I17) 
and the angular momenta of the noncondensate 

L iU) = e z ■ \ J dr[{u*(r x p)u q 

<m<kSt7<0'+ 1 ) 
+v q {r x p)v*}f(e q ) + v q (r x p)v*], (18) 

L 2 (e") =e,-\j ^(rxp){/( 6 (r )P )) 

to verify the mutual consistency of the Popov approxima- 
tion and that of the semiclassical approximation. Above, 
gi and L\ are obtained within the Popov approximation, 
while <?2 and L2 are obtained within the semiclassical ap- 
proximation. Results of these two approximations are 
consistent with each other for e > 5 ftwtr- 

The core-localized mode is mainly affected by the par- 
ticle densities inside the core, while g\ and <?2 account for 
the densities over the whole area of the system. Since our 
main interest is the core-localized mode, we employ cal- 
culations with reduced accuracy for n-i and gi for larger 
x,y (around \x\ > 1.2i?TF or \y\ > 1.2i?TF)- This affects 
the plot for gi in Fig.^Ja), but it has little effect on e' as 
shown in the dependence of e' on the cutoff energy, 

"cut- 
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FIG. 3: Horizontal axis indicates the displacement Ar of the 
vortex core from the center of the harmonic trap, (a) The 
vertical axis is the angular velocity ui Iot of the rotating co- 
ordinates for temperatures T = 0.0,0.1,0.3, and 0.5 T c . The 
difference in the frequencies is less than 2% between a point 
at T = and a corresponding one at T — 0.1 T c within 
Ar < 0.4. The dotted line is the result of the TF approx- 
imation at 7 - Q. Our results for T = 0.1 T c are 16 - 19% 
above those of the TF approximation, (b) The lowest ex- 
citation energy. All of the energies in the finite-temperature 
(T = 0.1, 0.3 and 0.5) systems are positive. The energy of the 
lowest core localized excitation only turns negative for T — 
when the normal densities m and 712 are neglected. 



for a small displacement Ar. This relation remains valid 
for < Ar < 0.5i?TF within the accuracy > e 1 > 
— 0.0047iwtr in the present system, see Fig.[3{b). 
Equations (fTSft . (|2U|) . and (ED lead to the relation 



IV. THE SIGN AND DIRECTION OF VORTEX 
PRECESSION 

The angular velocity aj prec of the precessional motion 
of a vortex is related with that of the rotating frame w ro t 
through 

Wproc + W rot = COnst. (20) 

within the GP equations. The core-localized excitation 
has the lowest excitation energy e' within the ran ge o f 
o; ro t we treat. Using variational Lagrangean analysis [12j, 
it can be shown that 

e' = at w prec = (21) 



e' = hq 6 Lu prcc . (22) 

Therefore, the direction of the precessional motion and 
the sign of the core-localized excitation correspond to 
each other at T = 0. 

If the coordinate transformation in Eq. Q2()[| were not 
valid, the relation Eq. (|22ll between the excitation energy 
e' and the angular velocity oj prec of the precessional mo- 
tion would not hold. The next section describes how this 
occurs at finite temperatures. 

An easier way to disprove Eq. (|22|l is as follows. Figure 
13 (b) displays e' for u picc = 0. It shows that Eq. lt2"T|) 
is no longer satisfied for T > 0. Therefore, it becomes 
impossible to satisfy Eq. (|22|l . The sign of the lowest ex- 
citation and the direction of the precessional motion are 
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FIG. 4: (a) Spectral distributions gi(e) obtained within the 
Popov approximation and 32(e) obtained within the semiclas- 
sical approximation, (b) Angular-momentum distributions L\ 
and 1/2- These two approximations yield consistent particle 
numbers and angular momenta for a wide range of energies. 
The temperature T/T c = 0.5 and the angular momentum 
A{<j>) = 0.863. 



no longer related within the finite-temperature (Popov) 
approximation. 



V. PROCESSING AND ROTATING FRAMES 

We have presented density profiles of an off-centered 
vortex within the rotating frame. This frame mimics a 
system of atoms in a rotating trap although the deforma- 
tion of the trapping potential is not included explicitly. 
The vortices are precessing if we observe them from the 
nonrotating frame. But there is a restriction in the ve- 
locity of the precession. 

Let us compare the velocities of the noncondensate and 
the angular velocity of the frame w ro t at w proc = 0. Figure 
|5] (a) displays the local velocities on the x axis for the 
system. The velocities 



V 12 (r) 
V(r) 



Re[Pi(r) + P 2 (r)] 
m{n\(r) +n 2 (r)}' 
Re[P (r) + P x (r) H 



P2(r)\ 



mn(r) 



(23) 
(24) 



sities 

Po(r) 
Pi(r) 



<t>*(r)p<t>(r), (25) 
{ u q( r )P u i( r ) + v q( r )P v q( r )} f( £ i) 

q (e q <E cut ) 

+v q (r)pv* q (r), (26) 

1 / e KF (r,p) 

2 \ e(r,p) 

e[e(r,p)-E cut }. (27) 



1 



The velocity closely follows that of the rotating frame 
outside of the condensate (\x\ > Rtf) in Fig. Efa). ft i s 
straightforward to show that V(r) reduces to 



V(r) = LO Tot (-y,x,0) 



(28) 



which are particle-current densities divided by the parti- 
cle densities; they are defined using the momentum den- 



for large \r\ where the density of the condensate and nor- 
mal component below the cutoff are negligible (no+ni <C 
n). This indicates adiabaticity between the rotating trap 
and the normal gas. Within the Popov approximation, 
the angular velocity of the precessional motion of a vor- 
tex is restricted to that of a normal gas and a confining 
trap. 

We have considered static systems with a-> prcc = in a 
rotating frame with the angular velocity o-> rot . These two 
w's may be transformed between each other using the 
simple relationship between the stationary and rotating 
frames of reference in Eq. (|20|) . valid at zero tempera- 
ture. But taking into account the normal component, 
this relation is only valid when the normal component is 
rotating at the angular velocity w ro t. Varying w rot will 
change the density profiles through Eqs. © and ljl3|l . 
Then the coordinate transformation Eq. (|20|) does not 
hold. Adiabaticity required within the Popov framework 
restricts the recognizing of c<j ro t as the angular velocity 
of precession a; prcc . This is another reason why Eq. 
does not hold at finite temperatures. 



VI. DISCUSSION 

It is confirmed that the sign of the lowest excitation 
energy e' is, in general, unrelated with the direction of 
the precessional motion of a vortex within the Popov ap- 
proximation. Within the Bogoliubov approximation, the 
excitation energy e' of core-localized mode and the pre- 
cession frequency a> pr0 c are proportional to each other as 
shown in Eq. (J35J. The derivation of Eq. (|2*2)l shows that 
the angular velocity (and the direction) uj picc of the pre- 
cessional motion arises from the coordinate transforma- 
tion Eq. (|20|) . which does not have any explicit relation 
with the core-localized excitation. Therefore, the core- 
localized excitation is responsible for the inward/outward 
motions of vortex and not explicitly related to the pre- 
cessional motion. 

We think that this nature of the core-localized exci- 
tation does not change even at a finite T. However, a 
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FIG. 5: (a) Velocities of all the atoms (solid line) and the 
normal component (dashed line) along the x axis (y — 0). 
The velocity follows that of the rotating frame (dotted line) 
outside the condensate (\x\ > -Rtf). Temperature T = 0.5T C 
and angular momentum A(4>) = 0.863. (b) Density profiles 
along the x axis for comparison with (a). The plots of no 
(dashed line) and m + n,2 (dotted line) are equivalent to those 
on the x axes in Figs. |3{a-b). 

system described within the Popov approximation can- 
not undergo the coordinate transformation. Within the 



Popov approximation, the velocity uj [ot in Fig. |2ta) may 
be recognized as the angular velocity of the precessional 
motion only when the condition of adiabaticity is obeyed. 
Therefore, the relation between the sign of the lowest ex- 
citation energy s' and the sign of w prcc , i.e., the direction 
of the precessional motion of a vortex is broken within 
the Popov approximation. 

We obtained the precession frequencies of a vortex 
at finite temperatures. The vortex-precession frequency 
manifests [Fig. Efa)] the same tendency as a function 
of displacement Ar as the zero-temperature TF stud- 
ies The precession frequency increases as the dis- 
placement Ar and the temperature increase. The fre- 
quencies differ only little between those for T = and 
T = 0.1 T c , despite the large changes in the excitation 
spectra. Although a pinning effect 0, 1^3, caused by 
the normal component makes a significant difference in 
the excitation energy e', it has little effect on the preces- 
sion frequencies between the systems at zero and finite 
temperatures. 
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